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This work stresses the importance of the choice of the set of reference functions in the Gener- 
alized Multichannel Quantum Defect Theory to analyze the location and the width of Feshbach 
resonance occurring in collisional cross-sections. This is illustrated on the photoassociation of cold 
rubidium atom pairs, which is also modeled using the Mapped Fourier Grid Hamiltonian method 
combined with an optical potential. The specificity of the present example lies in a high density of 
quasi-bound states (closed channel) interacting with a dissociation continuum (open channel). We 
' demonstrate that the optimization of the reference functions leads to quantum defects with a weak 

' energy dependence across the relevant energy threshold. The main result of our paper is that the 

, agreement between the both theoretical approaches is achieved only if optimized reference functions 

^SJ ' are used. 

^ ^ PACS numbers: 

<: 

lO ' I- INTRODUCTION 

(N : 

The idea of quantum defect in atomic physics has been introduced many years ago for Rydberg atoms, to describe 
' their similarity and difference with the hydrogen atom: the presence of the ionic core in the Rydberg atom induces a 
phase shift, i.e. a defect, on the wave function of the Rydberg electron which occasionally travels through the core, 
■ compared to the pure coulombic wave function which is called the reference wave function. It is well known that the 
^ I usual Rydberg law for hydrogen level energies En = —Roo/n^, where Roo is the Rydberg constant and n the principal 
^ , quantum number, is then modified according to Ene ~ —Roo/{n — He)'^, where /i^ is the quantum defect associated 
qh' to a Rydberg series (referred to as a channel) characterized by the orbital angular momentum i of the outer electron. 
I— I, The Rydberg law can be easily extended to Rydberg electronic states of a diatomic molecule, by defining a quantum 
defect fJ.(^\^a+ (R) which now depends on the projection A of £ on the molecular axis, on the internuclear distance R, 
" kJ [ and on the internal quantum state a+ of the molecular ionic core. In a Rydberg system, the wave function of the 
1 Rydberg electron inside the region defined by a range of strong interaction with the ionic core is independent, 
I apart from a normalization factor, from the generally weak binding energy Enii\a+ (R) of the electron. The remaining 
. nodal structure of the wave function outside this region is only determined by the Coulomb tail of the electron-core 
CO ■ potential. The quantum defect is basically independent of Eni\a+{R), so that it can be extrapolated through the 
, energy threshold to describe electronic continuum states, i.e. ionization. In most systems, several Rydberg series 
' are present in the same energy range, which most likely interact together. Such situations are described with the 
Q^ ] Multichannel Quantum Defect Theory (MQDT)[ll, 0, which involves quantum defects for each of the p channels 
^> • and coupling parameters between them, constituting the so-called Y-matrix with dimension p x p 

The quantum defect concept is very appealing in the context of ultracold gases, as it expresses the same idea than 
the scattering length for clastic collisions between ultracold atoms. Such a collision induces a phase shift on the 
wave function of the relative motion of the atoms, created by their strong interaction when they lie together within a 
^ . distance Rq. The scattering length is an effective parameter directly related to this phase shift, and characterizes the 
short-range interaction between the atoms. Just like in the previous case, the wave function inside i?o does not depend 
- apart from a normalization factor - on the relative (weak) energy of the atoms, while beyond i?o, it is controlled by the 
long-range interaction of the atoms, expressed as a multipolar expansion Ck/ R^ of the potential energy. Several 
authors worked at designing a Generalized Multichannel Quantum Defect Theory (GMQDT) to treat interaction 
potentials with non-coulombic asymptotic behavior [J, [1, 0, Q , in order to represent the vibrational spectrum of small 
molecules close to their dissociation limit, or the collisions between atoms for energies just above the dissociation 
threshold. In contrast with the pure Coulomb field problem (defining electronic Rydberg series), there are an infinite 
number of possibilities for the choice of the reference functions for the molecular problem (involving vibrational series). 
Therefore the central issue is to find reference functions adapted to an arbitrary long-range variation of interaction 
potentials, in order to rigourously define quantum defects and relevant channel c oup 

lings B 9- Recently, GMQDT 

has been used to study Feshbach resonances in ultracold atomic binary collision [l^, [ll| . These studies have shown 
that GMQDT accurately recovers close-coupling (CC) cross sections results and brings more physical insight on 
the underlying molecular processes. 

Using the same approach than in ref. pd| , we focus our attention on the characterization of the position and the width 
of resonances in the photoassociation of ultracold atoms, which strongly depends on the choice of reference functions. 
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We demonstrate that the optimization of the reference functions [l^, [3l based on a procedure introduced in 
ref.[l^, leads to quantum defects with a weak energy dependence across the relevant energy threshold, just like in 
the standard MQDT situation. We compare our results to numerical coupled-channel calculations performed in the 
framework of the Fourier Grid Hamiltonian method [l6j combined with an optical potential approach for the resonance 
characterization (l7j . The main result of our paper is that the agreement between the both theoretical approaches is 
achieved only if optimized reference functions are used in the MQDT framework. 

Our paper is organized as follows. We first present in section |TT] the chosen physical situation, i.e. the photoassoci- 
ation of ultracold rubidium atoms into mixed electronic molecular states coupled by spin-orbit interaction. We recall 
next the main features of its numerical description using the Mapped Fourier Grid Hamiltonian (MFGH) method with 
an optical potential (section IIII[) , yielding vibrational energies and predissociation line widths for high- lying vibra- 
tional levels of the coupled electronic states. In section [TVl we present a summary of the GMQDT treatment required 
to calculate the position and width of resonances. We compare in section|V]the results obtained by the two theoretical 
approaches for both the discrete and the resonant spectrum. Finally we discuss the validity of the GMQDT parame- 
ters defined in previous analysis of this system, in order to shed light on their predictive power for future experiments. 
Atomic units will be used for distances (ao=0. 052917720859 nm) and energies (2i?oo = 219474.63137054 cm"^), except 
otherwise stated. 



II. THE {A, h) COUPLED STATE SYSTEM IN THE RUBIDIUM DIMER 

The coupling mediated by spin-orbit interaction of the and Ir'Ii-u electronic states (hereafter referred to as 

the A and h states, or as the {A^ b) coupled system) correlated to the lowest excited ^5* P dissociation limit in 
alkali dimers, represents a well-known case for the breakdown of the Born-Oppenheimer approximation which strongly 
perturbs their spectroscopy. This is particularly true in heavy alkali dimers like Rb2 [iSl . Il9|. which required novel 
methods to deperturb the spectra recorded in high-resolution molecular spectroscopic studiesj20l [2l] | , where almost all 
rovibrational levels exhibit a mixed singlet/triplet character up to the dissociation limits. This pattern has also been 
observed in heteronuclear systems involving one heavy alkali atom like NaRb [13, |2^ , NaCs , and RbCs [2^ . In the 
context of the photoassociation (PA) of ultracold rubidium atoms, this strong interaction results into perturbations 
and predissociation resonances which have been observed for the first time in re f. 12611 . Such perturbations have been 
further studied exp erimentall y in the PA spectra of Rb2 [13, HI] , and Cs2 [1^, l30j |. and modeled by various means 
[28l . [29l . [sol . IsU . \3% m, [13, [13, [3^. In front of this intense activity, the {A, b) coupled system in Rb2 appears as a 
suitable test case for the present study. 

In Hund's case a representation, the relevant molecular potential curves Va{R) and Vb{R) of Rb2 are associated 
to the A^E+ and b^Uu electronic states correlated to the lowest excited dissociation limits 5 ^S" -f 5 ^P. Due to 
their spin-orbit coupling, they define a subspace of 0+ symmetry in Hund's case c according to the corresponding 
interaction Hamiltonian 1^8° (0+): 



vi:)(o^) 



H-A-(i?)/2 
VV^A-JJiR) 



^A°JJ{R) 
Va 



(1) 



At large distances, both i?-dependent spin-orbit coupling terms A°^{R) and A°JJ {R) reach the atomic value 
2A£'so/3, where AiS^o = 237.6 cm~^ is the fine-structure splitting of the rubidium 5p atomic level. In the fol- 
lowing, we neglect the atomic hyperfine structure, as well as the rotation of the molecule, and we restrict our study to 
the ^^Rb isotope. Furthermore, as we focus our attention on the comparison between two approaches, we will assume 
A°^{R) = All^ [R] ~ Aso{R) in the following. Equation ([T]) is rewritten in the asymptotic basis where the atomic 
spin-orbit interaction is diagonal (i.e. at i? — > 00): 



VA/i + 2Vb/i^A,o{R) y^{VA~Vb) 1 

y^{VA-Vb) 2VA/i + Vb/i + Aso{R)/2\ ^' 

Equation [2] is convenient, as the diagonal elements now converge towards the 5^5 -I- 5^Pi/2,3/2 dissociation limits 
(hereafter referred to as the P1/2 and ^3/2 asymptotes), and the off-diagonal coupling terms vanish at finite distances. 
These diagonal elements and their coupling are drawn in Figure [T^, and represent the diabatic representation of our 
problem, with two potential curves encountering a real crossing around lOao, coupled by an off-diagonal term reaching 
its maximum around the same value. The diagonalization of V°q(OJ) at every R yields two Oj potential curves in the 
adiabatic representation (Hund's case c), with an avoided crossing around lOa n as well (Figure [Hd). We used here the 
same data than in ref.jsj]. In standard photoassociation experiments [26l. [27l. l37j . bound levels located below one of 
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FIG. 1: (a) Diagonal matrix elements (full lines) and coupling matrix element (dashed line) of the interaction matrix of 
{diabatic representation) for the rubidium dimer; (b) Eigenvalues of the VJo(Oj) matrix {adiabatic representation). For 
illustration purpose, bound levels and predissociating resonances are pictured with full and dashed horizontal lines, respectively. 
The energy origin is taken at the 5 + 5 asymptote (horizontal dotted line). 

these asymptotes are populated. Here we are interested in eomputing energies for truly bound levels below the P1/2 
limit, and predissociated levels lying between the P1/2 fmd P3/2 limits. 

III. THE MFGH METHOD WITH OPTICAL POTENTIAL 

We briefly recall below the main features of the MFGH approach combined with an optical potential (MFGH+OP), 
which has been extensively described in several previous papers of the group [H, [13, S, S, [s^ . We want to solve 
the system of two coupled Schrodinger equations for the {A, b) system expressed in the diabatic representation, which 
is conveniently written in a matrix form according to: 



1 

2n dR^ 



I + i?I-V-(i?) 



^iE,R) = 



(3) 



where I is the identity matrix, E the total energy, and fi the reduced mass of the rubidium dimer. In the perspective 
of the next section, we express the matrix V"^ in eqOwith the compact form: 



VooiR) Voc{R) 
Vco{R) Vcc{R) 



(4) 



where the indexes o hold for open and c for closed, related to the P1/2 and P3/2 asymptotes, respectively. The 
solution ^{E,R) of eq. [3] is a two-component wave function 'if{E,R) = [ipoiR, E),ipc{R, E)) with radial wave 
functions ipo{R,E) and ipdRyE) associated to the open and closed channels, respectively. 
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In the Fourier Grid Hamiltonian (FGH) approach [l^, [3l| , we define a grid of length L with N equally spaced points 
in the R coordinate, associated with a set of N plane waves (jikiR) = exp{2TTikR/ L), k = —{N/2 — 1), ...,0, ...N/2 in 
the momentum space. Radial wave functions are then expressed as an expansion '^(i?) = X^iLi '4^{Ri)(t^i{R) ^ in which 
expansion coefficients are the values of the wave functions at every grid point. The potential energy operator V^^ is 
represented by a 2N x 2N matrix, composed of: (i) one N x N diagonal block per channel, each of them containing 
only diagonal elements equal to the potential energy of the channel at each grid point; two identical NxN ofF-diagonal 
blocks, themselves also diagonal and equal to the coupling term at each grid point. The kinetic energy operator is 
represented by a matrix T which is block-diagonal, with one dense block for each of the o and c channels. In order 
to save grid points when calculating bound levels close to the dissociation limits, i.e. with large spatial extension, or 
predissociating levels, we define a spatial grid with a variable step size determined by the function s{R) which maps 
the variation of the local classical kinetic energy of the radial motion in the coupled state (MFGH method [3^): 



^2^l [E^-V,nf{R) + E^] 

where Vinf{R) is the curve built from the lowest of the Voo{R) and Vcc{R) potential energy values at every R, Eo 
the energy of the lowest asymptote, i.e. the P1/2 one, and Eoo the range of dissociation energies above Eo to be 
explored in the numerical application. 

The diagonalization of the full hamiltonian matrix H = T + V^^ yields 27V eigenvalues and eigenfunctions for the 
coupled system, the latter having themselves 2N components in this representation. In the energy range between the 
two dissociation limits P1/2 and P3/2, levels from the upper (closed) channel interact with the dissociation continuum 
of the lower (open) channel, inducing predissociation. Following rcfs.|33.[38|. we add a purely imaginary potential (or 
optical potential) Vopt {R) to the diagonal Voo term. We chose the expression proposed in ref . [l^l : 



Noptexp 



2Lopt 

R ~ Ropt 



(6) 



where the recommended value for the normalization factor is Nopt =13.22. The distance Ropt characterizes the 
position of the optical potential, which has to be carefully chosen, depending on its amplitude Aopt and its range Lopt • 
Indeed, Vopt{R) must be placed at distances well outside the range of molecular potentials, in the region where the 
potential energies and the couplings arc negligible. We then ensure the convergence of the accumulated phase of the 
predissociatcd levels, and of their width. Values of Lopt = 40ao and ^opt=0. 00004 a.u. were found satisfactory, while 
we placed the optical potential at the edge of the grid, i.e. by varying its position Ropt between 40ao and 140ao, with 
an upper bound for the grid ending between 80ao and ISOoq. 

The diagonalization of the resulting complex Hamiltonian yields complex eigenvalues Ek = E^ — iTk/2^ with 
k = 0,...,2N. Below the P1/2 energy, all cigencnergics correspond to bound levels with pure real energies Ei,. 
Between the P1/2 and P3/2 asymptotes, the main task is to locate, among all complex E values, those which are 
associated to quasibound, or resonant levels with energy Er, and a finite predissociation lifetime l/F^. This can be 
done by several ways, illustrated in Figure [2] As in ref. [ssj . we rely on a stabilization procedure, which consists in 
increasing the size of the grid, and therefore of the density of discretized states representing the continuum. The 
imaginary part of resonant eigenvalues Er, i.e. their width F^, should converge with increasing grid size, in contrast 
to those associated to continuum states (see Figure [2ji). We see that resonances are identified over most of the P1/2 
and P3/2 interval, with L =80 a.u. and Ropt =40 a.u.. Increasing the size of the grid indeed provide converged width 
over the entire energy range. We note that the use of the mapping procedure is particularly appropriate under such 
a circumstance where large grids are involved. In Fig. [5)3, the width of the resonances is well defined for resonances 
with energy lower than the chosen Eoo value. We see that all resonances located between the two asymptotes are 
well described only if E^o matches the atomic spin-orbit splitting. In other words, the dissociation continuum has to 
be properly handled by the chosen grid step to obtain converged predissociation width. Finally, another convenient 
analysis of the results is provided by starting from the Breit-Wigner law for the phase shift: 



which allows defining the classical time delay [40|, |4l|, \^ : 



QE ^{E- Ek? + (rfe/2)2 ^ ' 
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FIG. 2: (a) Imaginary part Tk/2 of the solutions of eq|3] yielded by the MFGH+OP method, as a function of their real part 
Ek, between the P\/2 and P3/2 asymptotes. Einfty is chosen to match these limits, and Lopt ~ 40 a.u. and ylopt=0. 00004 a.u.. 
Full black circles: L =80 a.u. and Ropt —40 a.u.; Open blue triangles: L —160 a.u. and Ropt —120 a.u.; Open red circles: 
L =180 a.u. and Ropt =140 a.u.. Predissociated levels are assigned to the converged Tk/2 values, (b) Same as (a) with L =160 
a.u. and Ropt =120 a.u., for different values of the Eoo =0 (black open circles), 50 cm~^ (blue closed triangles), 238 cm~^ 
(red open diamonds), (c) Time delay which units???? associated to the complex eigenvalues, calculated according to eq|8] 
for Eoo =8 cm~^ (blue dashed lines), and Eoo =238 cm^^ (red full lines). Several resonances are identified with vertical red 
dotted lines, for illustration purpose. 



where the summation over k concerns complex eigenvalues with real part larger than the energy of the P1/2 
asymptote. From Figl^t, we clearly identify resonances by their larger time delay, compared to the continuum states. 
As expected, it is crucial to match Eoo to the energy interval between P1/2 and P3/2 to describe time delays for all 
resonances. 
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IV. GMQDT TREATMENT 

We focus on a system of one open channel and one one closed channel (labeled by o and c indexes respectively) 
coupled together by a coupling vanishing beyond R > Rq. Following the previous section, they correspond to the 
electronic states correlated to P1/2 {open) and P3/2 {closed). The Schrodinger equation has two linearly independent 
solutions, which can be written beyond Rq according to the following matrix form: 

^{E, R) = f{E, R) - g{E, R)Y{E), R > Ro (9) 

In this equation, all boldfaced symbols are squared matrices: '4'(£', R) is the set of two independent solutions, and 
f and g are diagonal matrices. Their respective diagonal elements {io{R, E) , ic{R, E)) and {go{R, E) , gc{R, E)) are 
the sets regular and irregular energy- normalized reference functions associated to the open and closed channels. The 
matrix elements of Y(_E) are the phase shifts induced by the short-range interaction, parameterized by the quantum 
defect matrix according to Y(£') = tan/x(-E'). This equation clearly emphasizes the link between the choice of 
reference functions, which is not unique, and the definition of the quantum defect. For instance, when the interchannel 
interaction is weak one obvious choice are the solutions of the Schrodinger equation for the individual open and closed 
channel, involving the Voo and Vcc diagonal matrix elements of V"^. In contrast, when the interchannel interactions 
is strong a better choice is the functions associated with the adiabatic potentials resulting from the diagonalization 
of V°*. Other possibilities have been discussed in the literature, like the analytic reference functions specific to or 
long-range potentials [i^, [3] , or the numerical Milne solution [4^, [4^ of the long-range potential as proposed by 
Greene and coworkers d]. The common property of all these reference functions is that the associated Y{E) matrix 
of equation ^ does not present any discontinuity across the dissociation threshold. Whatever the choice of reference 
functions is, the method yields the bound states energies and cross-sections, while the characterization of resonances 
strongly depends on this choice, as we will see below. 

Let us start with the GMQDT treatment described in ref. [i^l, for which we outline its main steps below. Assuming 
a strong interaction between the open and closed channels pictured in Figure [TJi, we determine the reference func- 
tions f and g"'**'* associated to the individual adiabatic channels of Figure [T}d with the Milne phase-amplitude 



method |45j . In order to avoid the tricky numerical evaluation of the non-Born-Oppenheimer coupling between the 
adiabatic channels, the set of coupled equations of eql3]) in the diabatic representation is solved by propagating the 
related diabatic logarithmic derivative matrix L'''" of the wave function ^ using the renormalized Numcrov method 
implemented by Johnson If M(i?) is the eigenvectors matrix of Hgo, it can be shown that the corresponding 

adiabatic logarithmic derivative matrix L"'*'" is given by: 

L°*°(i?) = M{R) L*"(i?) M*(i?) (10) 

The coupling between the adiabatic channels is localized around their avoided crossing (see Figure [T)d) . Outside 
the interaction region, i.c for R > Rq = 13a. u. in the present case, the Y'"''" matrix elements no longer depend on 
R. Therefore equation ([9]) holds, and Y"''*" is extracted according to: 



Y 



adia 



L"*°(i?o)g"*°(i?o) - (g''*°)'(i?o) L"*°(i?o)f"*''(i?o) - (f"*°) (i?o) (11) 



-1 



adia f \radia/T^ \ fradia\ i 



The extracted Y'^''*" matrix involves both the closed and open channels, and relevant asymptotic conditions are 
applied according to the spectral range of interest. For energies below the P1/2 asymptote, bound states energies E^ 
are determined by solving the equation: 



\t^nu,{Eb)+YT''{Eb)\^Q (12) 

where the vertical bars hold for the determinant of the matrix. In this expression, tan^^c is the closed block of the 
diagonal matrix associated to the accumulated phase i'{E) of the adiabatic channels, which is determined numerically 
using the Milne method. 

In the energy range between the two asymptotes, the long-range behavior of the closed and the open channels 
are different. As implemented in ref.[l|, the scattering matrix S is obtained by applying the relevant asymptotic 
conditions to first extract the open block of the reaction matrix K defined by: 

'Koo{E) = Ylf" - Y^jf " [tanv,{E) + Y'^^f''] Y^^^'" (13) 
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FIG. 3; Energy variation of the Yq"'*"' (crosses), Y'^c'*"' (plus signs), and Fq"/"" (open circles) matrix elements in the region 
around the P1/2 threshold. The dashed lines illustrate their linear extrapolation from energies below P1/2 across the threshold. 
The off-diagonal element of the optimized matrix Y°^^ (full circles) presents a weaker energy dependence. 



The S matrix is then restricted by definition to the open channels: 



where {^{E))o, are the usual shift of the reference functions with respect to the Bessel reference functions, for 
the open channels. We have omitted the energy dependence of Y matrix on purpose, to emphasize that the energy 
dependence of the scattering matrix is mainly governed by the energy variation of the parameters ^{E) and v{E) 
and that the resonances arise from the pole structure of equations (|13ll4p . The Y matrix elements are expected to 
slowly vary with the energy, provided that the choice of reference functions is appropriate. On the numerical side, 
they must be evaluated on a thin energy mesh while the Y matrix elements only require to be evaluated on a coarse 
mesh followed by a spline interpolation. 



In the present model, we consider the two-channel problem described in section |TT1 so that there are only one open 
channel and one closed channel. Most of the previous GMQDT equations simplify, as the open and closed blocks 
reduce to scalar quantities. First, the energy variation of the Y'"*"* matrix elements is displayed in Fig[3]in the region 
around the P1/2 threshold. With the present choice of reference functions associated to the adiabatic representation, it 
is clear that their extrapolation from the discrete range (below the threshold) through the energy continuum between 
the dissociation limits is rather hazardous. 

Before addressing the issue of the choice of the reference functions, it is useful to analyze the results obtained with 
the GMQDT approach in the adiabatic representation, compared to those of the MFGH method of Section IIIII The 
bound state energies Et of eq[T2]are reported in Figure Sl^b) , compared to those obtained by MFGH over a 80 cm~^ 
energy range below the threshold. In this figure, each energy value is given an ordinate corresponding to the 
squared modulus of the wave function for the P1/2 channel, hereafter referred to as the weight of the lowest channel 
in the total two-channel wave function. The influence of the bound levels of the upper channel due to the short-range 
coupling is clearly visible, inducing the minima in the weights of the P1/2 components. As expected, the agreement 



Soo{E) = exp(+i ^{E))o [lo + i Koo{E)] [!„ - i Koo{E)] ' exp(+i ^{E)) 



O 



(14) 



V. COMPARISON OF MFGH AND GMQDT RESULTS 
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FIG. 4: Weight of the lowest channel in the total wave function (see definition in the text) for bound states located below the 
Pi/2 limit, as a function of their energy position Et: (a) by the MFGH approach and (b) by the GMQDT approach. 



on bound state energies is perfect between the two approaches. In contrast, it is not surprising that the weights are 
not identical in the two representations, as the chosen channels are not defined in the same manner, as it can be seen 
by comparing the two panels of Figure [1] 

Next, the results in the energy range between the two dissociation thresholds where the resonances are expected 
to take place are interpreted through the S matrix elements, of more precisely through its eigenphases 5s defined 
as S = exp(2i 5s)- We reported in Figure [5] the quantity sin^ 5s (sec the full lines in the figure), which shows that 
its energy variation is so irregular that it is impossible to recognize any resonant pattern. This suggests that the 
dynamics of both channels is strongly coupled to each other. One can isolate the contribution of the resonances to the 
total phase shift Jg, by first setting the notation tan 5k = so that 5s can be extracted from eq. [TD 5s ~ £,o + 5k- 

If we neglect the channel interaction, the Koo reduces to Y^g*" and we can define a so-called background phase 
shift according to tanJ^g = Y'^f'"- . The resonant part 5r of the phase shift can be calculated as 5r = 5k — 5bg- The 
energy variation of sin 5r between the two asymptotes is displayed in Figl6K, where the regular series of resonances 
can now be identified from the energy of the P1/2 limit until about 100 cm~^ upwards. One can note that even in the 
low energy part the resonance profiles do not correspond to a Lorentzian shape, as we would expect from an isolated 
resonant contribution. A gradual increase of the profile asymmetry with increasing energy is visible, decreasing again 
close to the P3/2 limit. For comparison, one can also extract this resonant contribution from the MFGH+OP approach 
(FigE^). Resonant energies Er and width are deduced as outlined in Section Hill which assumes that the related 
resonances are described by a Brcit-Wigner profile (eql?]) which can be expressed in an equivalent way as: 



-'■■'^-'^'- i^j,y!';r,/2)^ <'^) 

We see that in Figure [6] both approaches yield identical resonance positions up to about 160 cm^^ above P1/2. 
However, due to the asymmetric resonance profile, the comparison of width is not obvious. Actually, the resonance 
width can be directly extracted from the Y"'^*" matrix, using a procedure proposed by Lecomte [I3. The starting 
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FIG. 5: Energy variation of the total phase shift Ss between the P1/2 hmit (taken as the origin of energies) and the P3/2 
limit. Full lines correspond to the GMQDT calculation performed with the adiabatic set of references functions. Full circles 
correspond to the results obtained with the rotated set of reference functions. 
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FIG. 6: Energy variation of the resonant part of the phase shift 5r between the P1/2 limit (taken as the origin of energies) and 
the P3/2 hmit. (a) MFGH calculations; (b) GMQDT calculations. 
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point is to eliminate the open channels to build the block of the K matrix associated to the closed channel by applying 
Siegert aymptotic conditions, instead of doing the reverse as in eall3l (we omitted the adia index for simplicity): 

Kcc = Ycc - YcoYoo [I + Y^J Yco + tYco [l + Y^^] (16) 

Following the framework of configuration interaction as described by Fano [i^ , one can interpret all terms of this 
expression. The imaginary part is due to the interaction between the closed and open channel Y^o which give rises to 
resonance in the collisional cross-sections, and is proportional to the width of the resonance. The real part consists in 
two terms. The leading term is the direct coupling term Ycc between the closed channels. It is modified by the next 
term which describes the indirect coupling among the closed channels through their interaction with the open ones. 
This latter term represents the shift of the resonance relative to its zeroth-order position determined by the direct 
interaction term Ycc. The resonance energies Er are located by solving the determinental equation: 

tan z/c(Sr) + Y,, - Y,oYoo [l + Y^] Y,, ^ (17) 

which reduces to a scalar equation in our case with one single closed channel and one single open channel. The 
zeroth-order energy position of the resonance is obtained by canceling all interaction terms in eg 1 171 such that: 

t&ny,{E^°'^) + Y,, = (18) 

Using a first-order development of tani/c(£'r around E'r*'"' in eq |171 the resonant energy can easily be extracted as 
— -Bq — A = 0, with the energy shift: 

A = N-'\Eo)Y,o [{Yoo{I + Yl)-^] Foe (19) 
where N'^{E) is the usual QDT bound state normalization factor [l|: 

and is the usual QDT bound state normalization factor (Seaton 83). The associated width is then given by: 

= 2N-^{Er) Fco [(/ + Y^)] Yoc (21) 

The variation of the resonance widths with their energy Er is presented in Figure [71 as computed by both 
GMQDT and MFGH approaches. One can observe a good agreement between the two methods in the low-energy 
domain, while big discrepancies occur around 175 cm~^ above Pi/2- In this region, the resonance can no longer be 
considered as being isolated as their width become larger than the energy separation between consecutive resonances. 
Therefore the expression of eq |21l no longer hold, and the apparent narrow widths computed with GMQDT no longer 
reflect the actual strong coupling between the two channels For instance, their apparent widths in fig l6] are 
much smaller than the values given by Eql2T| which usually referred to as stabilization effect. The asymmetric 
(non-Lorentzian) profiles result from interferences between neighboring overlapping resonances. 

VI. OPTIMIZATION OF GMQDT REFERENCE FUNCTIONS 

The apparent failure of GMQDT to properly characterize resonances in the dissociation continuum actually arises 
from an improper choice of the set of reference functions. Indeed, we have seen in Figl3] that the intuitive choice of 
reference functions associated to the adiabatic channels induce large Y matrix elements with strong energy variation 
across the dissociation threshold. In other words, the non-adiabatic couplings are so strong that the adiabatic 
channels do not represent a good zeroth-order approximation in the present case. Several authors have shown that 
an appropriate rotation of the reference functions can improve such a situation. The basic idea is to define a new 
Y matrix where only the off-diagonal block is not zero. In such a representation no energy shift would occur on the 
resonance, i.e. the coupling between the open and closed channel would be minimized. Following ref . (isj . this goal is 
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FIG. 7: Energy variation of the resonance width given by Ea l21l as a function of the energy spacing of the resonant levels. 
Circles; MFGH; Plus signs: GMQDT. 

achieved in two steps. Two successive rotations of the initial adiabatic fc and gc (resp. /o and go) reference functions 
with a characteristic angle Oc (resp 6'o)are applied according to: 

= fc cos 9c - gc sin 9c (22) 
gc = fcsm9c + gcCos9c (23) 

with a similar equation for /o and go with the angle 9o. After each rotation, the new matrix Y"""* is given by: 



-^rot 

oo 
~\rrot 



CO cc 



-\j-rot 

oc 
~s^rot 



Y Y 

^ oo ^ oc 

Y Y 

^ CO ^ cc 



-co ^ cc 



COS 9o 

cos 9r 



sm6 




smb 






sin 6 



cost 






sin 6 


cos 9r 



(24) 
(25) 



with 9o = for the first rotation, and 9c ^ for the second one. First, the energy-dependent angle 9c is chosen 
such that the real part of the new K™* matrix associated to Y''°* vanishes, leading to: 



tan 29c = 



2Re{Kcc) 



1 - /m(K,,)2 - Re{Kcc)^ 
In this rotated reference set, eq[T7]used to locate the resonance simply reduces to: 



(26) 



taniiydE) + 9c) = 

For the second step, Y°p* is set to zero, which is fulfilled for the value 9o defined by: 



(27) 



tan 9o = Y, 



rot 
oo 



(28) 



It can be easily derived from Eql25] that this also ensures that Y°p* is zero. Hence, the Y"^* matrix is well 
transformed into a matrix with non- vanishing off-diagonal blocks. 

The energy variation of this single matrix element Y°p* is reported in FiglS] One clearly observes that the magnitude 
of the matrix element, associated with the rotated reference functions, is now smaller than unity. Moreover its variation 
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FIG. 8: Energy variation of the resonant part of the phase shift 5r between the P1/2 limit (taken as the origin of energies) and 
the P3/2 limit, (a) MFGH calculations; (b) GMQDT calculation with optimized reference functions. 

with the energy is much weaker than the variation observed for the ofF-diagonal element associated with the non- 
rotated reference functions. 

The energy variation of eigenphases of the scattering matrix determined with the optimized reference set is presented 
in FigO above. It must be noted that this result is obtained with eqs[T3]in which must be replaced by + ^o, 
Uc by i^c + and the Y matrix by its optimized counterpart. At first glance, as it should be, this variation exactly 
matches the one obtained with the non-rotated reference set. Within this rotated frame, the analysis of the resonant 
contribution 6r to the scattering matrix eigenphase 63 reported in Figure [51 shows that the resonances now appear 
with the expected Lorentzian profile. Moreover both theoretical approaches leads to the same results for the resonance 
location and width. 

A more quantitative comparison of the width of the resonances is reported in Figl^l The agreement between 
GMQDT and MFGH shows that the Fermi golden rule is recovered in the rotated set of reference functions. Al- 
though the collision cross-section energy variation is independent of the choice of reference functions, this application 
demonstrates how this choice is important for the analysis of resonances arising in the cross-section.. 

VII. CONCLUSION 

In this paper we emphasized on the importance of the choice of reference functions in the GMQDT treatment of 
collisional resonances, illustrated by the case of the photoassociation of cold atoms. In contrast with the case of 
interacting Rydberg series, where the quantum defects are always associated with the regular and irregular Coulomb 
functions, the present two-channel molecular problem allows several choices for the reference functions. We demon- 
strated that the widths of the resonances can be properly assigned provided that the reference functions are chosen in 
a way which minimizes the interaction with the related channels, and which ensures a weak variation of the associated 
quantum defect with the energy. Therefore the width of the resonances can be safely deduced after the extrapolation 
through the dissociation threshold of the generalized quantum defects computed for the bound state spectrum. We 
want to stress that both set of reference fimctions lead to the same collisional cross-sections because they correspond 
to a different balance between the background (governed by the reference functions) and the resonant contributions 
to the S matrix eigenphase. They only lead to a different analysis of the resonance showing up in the cross-sections. 
The non-rotated set corresponds to a strong coupling situation for which it is difficult to have predictable physical 
insight, while the rotated set corresponds to a more familiar and predictable weak coupling situation. 
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FIG. 9: Energy variation of the resonance width obtained from MFGH+OP (open circles), and from GMQDT using rotated 
reference functions (plus signs), as functions of the energy spacing between the resonant levels. 



With this in mind, it is worthwhile to revisit the interpretation of the experimental results reported in ref.[28[ 
concerning the photoassociation of cold ^^Rb atoms into the 0„ + {A, b) coupled states. In Figure \T0\ we reported the 
positions of the recorded PA lines below the P1/2 asymptote, and the energies of the corresponding levels computed 
with the GMQDT approach using the rotated reference functions above. We see that the molecular potentials used for 
the calculations provide energies in reasonable agreement with the experimental values, even if there is room for their 
improvement. More specifically, within the Lu-Fano approach (50| used in refs.[28l [23|. the underlying assumption 
is that the reference functions are associated to the long-range adiabatic 0+ potentials varying as even if these 
functions are not explicitly derived. This agreement simply means that our potential curves exhibit a long-range 
behavior which is close to the actual one. In this figure, the computed weight of the wave functions on the P3/2 
channel is used to yield an intensity to the lines. The resonance induced by the upper channel corresponds to the 
highest value of this weight. Over the displayed energy range, two of such resonances appear, separated by about 
7 cm~^. We note however that for each of these resonances, we predict the existence of one line undetected in the 
experiment, marked with a star in Figure fTOl They correspond to pairs of one P3/2 and one P1/2 level, which are 
the more strongly coupled to each other. It seems that only one partner of such pairs has been identified in the 
experiment. 

Such a representation is helpful to define the energy range of the P1/2 bound spectrum which is perturbed by the 
interloper belonging to the P3/2 channel. Indeed, the double-sided arrows in the figure suggest that the half- width of 
this perturbation is equal to about 0.5 cm~^, in agreement with the computed half- width of the very actual resonances 
lying just above the P1/2 threshold reported in Figure [9l This value is also in agreement with the energy width of the 
variation of the quantum defect around the resonances given in ref. [2^. In this latter paper, the authors extrapolate 
their model through the dissociation limit, which is indeed justified by the present analysis. For the first resonance 
above P1/2, their quantum defect indeed varies over a half-width of about 0.4 cm~^. However at the predicted energy 
position in their signal, they assigned a broader (2 cm~^) feature to this resonance. Therefore we suggest that this 
feature could be revisited experimentally, as well as higher resonances, hopefully to rule out any unwanted contribution 
to its broadening. 

Finally, predissociation resonances over about 35 cm^^ below the P3/2 threshold of ^^Rb2 have been observed in 
ref . [27t . Their spectrum show a series of resonances with decreasing half- widths, which can be estimated from their 
figure as varying from about 0.8 cm~^ down to about 0.02 cm~^ close to P3/2. This trend is in agreement with the 
one visible in Figure [HI where we predict a quasi-linear decrease of the half-width in this region, from 0.06 cm~^ 
down to almost zero. A more accurate determination of the half- width from raw data of ref. [131 is desirable to better 
compare its amplitude and variation with our prediction. It would represent a test of the spin-orbit coupling used in 
the calculation. In this respect, a new perspective for the present analysis is offered by the recent determination of 
the potential curves and spin-orbit coupling of the {A, b) coupled system in Rb2 , extracted from the compilation of 
several high-resolution spectroscopic data [l^. Indeed, the present calculations could be run again with these new 
accurate data, in order to predict the predissociation width over the whole energy range between P1/2 and P3/2, and 
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FIG. 10: Energies of the bound states of the *^Rb2 0j(y4, fe) coupled states computed within the GMQDT framework with 
optimized reference functions (see text), compared to experimental results of ref . [28] . below the P\/2 asymptote taken as the 
origin of energies. Intensities of the vertical lines represent the weight of the P3/2 channel in the GMQDT results. Experimental 
results are weighted by the same quantity, for convenience. Grosses indicate missing data in the experimental signal due to 
technical reasons. Stars indicate levels predicted by our model, which are not detected in the experiment. Double-sided arrows 
indicate the approximate energy width of the perturbation induced by levels of the P3/2 channel in the P1/2 channel. 



hopefully confirm the interpretation of the PA results of ref. [27 
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